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Abstract 

We study the behavior of energy levels in two dimensions for exotic atoms, i.e., when 
a long-range attractive potential is supplemented by a short-range interaction, and compare 
the results with these of the one- and three-dimensional cases. The energy shifts are well 
reproduced by a scattering length formula 6E = A^/ ln{a/R), where a is the scattering 
length in the short-range potential, Aq/{2 7:) the square of the wave function at the origin 
in the external potential, and R is related to the derivative with respect to the energy of the 
solution that is regular at large distances. 

1 Introduction 

Hadronic atoms give valuable information about strong interactions at low energy. For a review, 
see, e.g., [1|. They have also motivated several studies on the behavior of the energy levels in 
a Schrodinger operator, with a potential Vi + \ V2, where Vi dominates at large distances, but 
is superseded by V2 at short distances. The case of exotic atoms corresponds to a world with 
three dimensions, where Vi = — 1/r (as a negatively-charged hadron orbits near the nucleus 
and is almost unscreened by the remaining electrons, if any), and V2 describes the short-range 
hadronic interaction. But the situation is far more general, and many features do not depend 
on the Coulomb character of Vi. Nevertheless, we shall use the word "exotic atom" for such a 
system, "atomic" for the energy domain of the eigenstates of Vi alone, and "nuclear" for any 
typical energy within V2 alone, for the sake of simplicity. 
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The spectral problem of exotic atoms [[T]-[3j differs significantly from the ordinary perturba- 
tion theory, for which an expansion of the eigenenergies in powers of A is attempted. For exotic 
atoms, the energies for Vi + \ V2 are often very close to the ones for Vi alone, but perturba- 
tion theory usually does not hold. For instance, if A V2 is an infinite hard core of small radius, 
the energies are slightly shifted upwards, but the ordinary perturbative expansion diverges al- 
ready at the first order. The proper concept here is "radius perturbation theory", as described by 
Mandelszweig [|4|. 

In this paper, we discuss how exotic atoms behave in = 2 dimensions. It may be noted 
that the study of exotic atoms in d = 1 is more straightforward, and already discussed in the 
literature [[3j. The d = 2 case is more delicate. The leading order term for the energy shift is 
easily identified, and linked to In a, where a is the scattering length in the short-range potential. 
As in the d = 3 case, the overall coefficient is the square of the wave-function at the origin in the 
external potential. However, the scale regularizing this leading term, i.e., the radius R leading 
to In a — )■ \n(a/R) is not immediate, but it can be derived from a matching of the solution of Vi 
which is normalizable to the asymptotic solution emerging from the short-range term A V2. 

The case of d = 2 dimensions is rather special in spectral problems, as it corresponds to 
the largest value of d for which an attractive potential, however weak, always holds at least one 
bound state, see, e.g., [[5}[6|. [^Hence, for < 2, if V2 is attractive, AV2 immediately develops its 
own bound state, which becomes the ground state of the Hamiltonian. However, this process is 
less effective for = 2 than for = 1, and the spectrum, as a function of A evolves more slowly. 

This paper is organized as follows. In Sec. [2] and Sec. |3| we give a brief reminder about 
the cases of = 3 and d = 1 space dimensions respectively, with particular emphasis on the 
phenomenon of level rearrangement and on the scattering length (hereafter referred to as SL) 
formula for the energy shifts. In Sec. |4| we present the results for the case of = 2 dimensions. 
The theoretical framework is presented in Sec. |5j before the final discussion in Sec. |6j 

2 Exotic atoms in three dimensions 

There is an abundant literature about exotic atoms in three dimensions, motivated by experi- 
ments with pionic, kaonic and antiprotonic atoms [[T||2|. The simplest model consists of a two- 
component potential 

V1 + XV2, (1) 

where Vi is a long-range interaction with one or several bound states. Genuine exotic atoms 
correspond to Vi{r) oc — 1/r. The second term, with an explicit strength A introduced for the 
ease of the discussion, accounts for the short-range interaction. The main results are: 

• the shift is usually rather small, although A V2 can be very large at short distances, 

• the shift is usually well described by the approximate formula 

(5E = E(A)-E(0) ~47r 10(0)1^0, (2) 
'More precisely, what is sufficient is that the integral of the potential over the whole space is positive. 
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where is the normalized wave function for A = 0, and a the scattering length in V2 alone. 
In case Vi is Coulombic, one recovers the well-known SL formula by Deser, Golberger, 
Bauman and Thirring, and Trueman [171 



En - Ei'^ _4a 

where B is the Bohr radius and n the principal number for the energy En^ in Vi alone or 
En in the total potential. Many improvements and further corrections to this formula have 
been discussed in the literature [|9l[T0l. 

When A is varied, the shift usually varies very slowly, except near the specific values Ai, A2, 
. . . , where the energy levels change very rapidly, and a level rearrangement occurs: near 
A = A„, the n^^ energy drops toward very large (negative) values in the nuclear domain, 
and is replaced in the upper part of the spectrum by the next level, which in turn is replaced 
by the next one, etc. An example is given is Fig.[T] Further examples are provided, e.g., 
in [[3|. The critical values A„ correspond to the coupling thresholds for which the short- 
range interaction A V2 starts supporting a first or an additional bound state. 




Figure 1: Level rearrangement in three dimensions. A wide and weak external square well 
V = — Ai 9(ri — r), with Ai = 4 and ri = 1, is supplemented by a short-range square well of 
increasing depth, A V2 = —A 6(r — r2), with r2 = 0.1, in units where h^/ {2fi) = 1, /i being the 
reduced mass. The first few energy levels are shown against A r"^. Bottom: first three levels. Top: 
magnification for ground-state alone, with the exact value (thick line) compared to the first-order 
perturbation theory and SL formula. 



3 Results in one dimension 

An example of spectrum of exotic atom in d = 1 is shown in Fig. [2} It consists again of a 
superposition of two square wells, the strength of the short-range one being varied. The main 
differences, as compared to the more familiar d = 3 case are: 

• As soon as A slightly departs from zero, the atomic ground state energy immediately drops 
towards the range of the nuclear energies. 
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As a coupling threshold in A V2 is reached and A further increases, a plateau is observed; the 
corresponding energy drops, and, by rearrangement, a upper level makes another plateau 
near the same value. This plateau in the sector of the even parity states, corresponds to an 
unperturbed energy level in the odd sector of Vi. Indeed, the orthogonality with the ground 
state forces a zero in the wave function near a; = 0, and mimics an odd state. 

The Deser-Trueman formula, if translated iox d = 1, reads 



5E ~ -2 



(4) 



The presence of the scattering length a in the denominator can be understood by dimen- 
sional analysis. Also, weaker the short-range interaction A V2, more flat the zero-energy 
wave function, and thus larger the scattering length a, defined (as for c? = 3) as the abscissa 
where the asymptotic zero-energy wave function vanishes. 
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Figure 2: Level rearrangement in one dimension. Same interaction and same notation as Fig.[T] 



4 Results for two dimensions 



The calculation can be repeated for the isotropic (i.e., azimuthal quantum number m = 0) states 
with d = 2. If the atomic spectrum is examined for increasing values of the strength of the 
short-range interaction, a pattern of level rearrangement is clearly identified, see Fig. [3} 

The behavior of the ground state is displayed again in Fig. |4} where it is compared to the 
d = 1 and = 3 cases. The trend is clearly intermediate between the plateau of (i = 3 and the 
immediate fall-off oi d = 1. 

For small values of A, we can easily identify the following behavior for the d = 2 energy 
shift 5E 



6E 



42 
^0 



(5) 



ln{a/R) ■ 

If one plots, as in the example shown in Fig. [5} —1/5 as a function of In a, one hardly distin 
guishes the exact values from the results of a linear fit. 
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Figure 3: Level rearrangement in two dimensions, for the first four levels. Same interaction and 
same notation as Fig. [T] 




Figure 4: Ground- state energy for = 1 (left), d 
the same interaction as in the previous figures. 



2 (middle) and 3 (right) dimensions, with 



As discussed below, ~ 27r |0(O) p and a is the d = 2 scattering length, as recently revisited 



[11 12 1 . The value of R is found of the order of magnitude of the "Bohr radius" of the wave 
function in the external potential, that is to say, the average radius. Its expression is derived in 
the next section. 




In a 



Figure 5: 1/ 5E against In a for the ground-state energy in the double square- well of Fig. |4| The 
linear fit cannot be distinguished from the exact results. 
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5 Derivation of the energy shift 



5.1 General formula 

There are many approaches to the SL formula for d = 3, and various corrections and generaliza- 
tions, see, e.g., [[3}|4||7 - 10 1 and references there. For the d = 2 case, the following simple minded 
derivation is just based on the matching condition between solutions of the Schrodinger equation 
that are regular at short and at large distances. 

For the sake of clarity, one can identify several approximations that are made when solving 
the bound-state problem in the potential Vi + A V2: 

1 . Vi dominates at large distances 

2. the energy E in Vi alone is a smooth function of the boundary condition enforced at r = 0, 

3. Vi and the energy term E u can be neglected at very small distances, where A V2 dominates. 
Let us start with Schrodinger equation for the external potential alone, i.e., 

-u" - ^ + Viu- Eu = , (6) 

where u is the reduced radial wave function and we are working in the units with h = 2m = 1. 
We denote h(E,r) the solution that is regular at infinity, i.e., h{E,r) oc ^/rK^ikr) at large 
distance, with E = — /c^ and Kq the usual Bessel function. The case of a confining interaction 
Vi is treated later. At short distance, this solution behaves as 



h{E,r) = B{E) y/r \nr + A{E) y/r + ■ ■ ■ , 
h{Eo,r)=AoV^+--- 



for the modified energy E and the unperturbed one Eq. The unperturbed energy corresponds to 
B{Eq) = 0, i.e., a solution that is regular as r — i- 0, and is normalized, leading to a real value 
Aq = A(Eo) at energy Eq, that can be chosen to be positive. For E ^ Eo'm the neighborhood of 
Eq, we impose that the solution remains normalized, i.e.. 



+00 

h{E,ryr = l. (8) 







By combining (|6]) for h(E, r) and h{Eo, r), one obtains the exact relation 



+00 



Ao B{E) = {E- Eo) / h{E, r) h{Eo, r) r (9) 

which gives for the energy shift SE = E — Eq a first relation 6E ~ i? Aq. It is rather precise. 
Indeed, if the solution is kept to be normalized as per ([8]), and if h(E, r) — )■ h{EQ, r) as E ^ Eq, 
the integral of h{E, r) h{Eo, r) entering ([9j) is also equal to 1, up to second order in 5E. 
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Now, A = Aq + Aq 5E + ■ • ■ , one can identify the short-range behavior of h{E^ r) with 
[In r — In a], to obtain 

B Ao + Aq 5E 



1 — In a 

which when combined with B c:^ 5E/Aq gives 



(10) 



/12 

5Ec^ 5 — ^. (11) 

In a + Aq 

where the denominator can be cast as In a — In R. This relation gives explicitly the link between 
the energy and the boundary condition at r = 0, expressed by In a, where a is the scattering 
length in the short-range interaction alone, in terms of the quantities Aq and Aq linked to the 
value of the unperturbed solution at the origin. 

The effective range correction to the scattering length approximation can be worked out ex- 
plicitly, but turns out to be very small in most cases. For a positive energy E = k'^, the m = 
solution to the scattering problem in A V2 alone is pT[ - [T3| 

U2{k, r) = [cot 6{k) Jo{k r) + Yo{k r)] , (12) 

where Jq and Yq are Bessel functions, and 

cot 6{k) = - [ln(a k/2) + -f] + I ro k^ + ■ ■ ■ (13) 

TT 2 

involving the scattering length a and effective range tq. This expression is easily translated for 
negative energies E = —k"^, and, anyhow, the range of E which is explored is small as compared 
to the typical energy scale in the short-range potential, and thus the solution coming out is safely 
approximated by its small E hmit 

M2 = ^A' (Inr + Ina =F TT ro -E/4) . (14) 

This means one can replace In a by In a + nroEo/A to probe the contribution of the effective 
range, which turns out negligible, provided the energy shift SE remains small as compared to 
^0. 



5.2 First example: double delta-shell 

To illustrate ( fTT] ), we consider as long range interaction an attractive delta-shell of strength gi 
and radius that can be set to i?i = 1 to fix the length scale. The solution can be worked out 
analytically, in particular ^ KQ{kr) is regular at large r and -Jr I^ikr) at small r. The delta- 
shell interaction imposes the continuity of the radial solution u near r = Ri and the proper step 
in its derivatives, to fix the unperturbed energy Eq. A second delta- shell can be implemented 
at r = R2 <ti Ri, leading to an explicit transcendental equation for the exact energy E, and 
shift 5E = E ~ Eq, to he compared to the simple approximate value 5E' = —A^/ Ina and the 
improved 6E" given by ( [TT] ). For gi = 1/2, R2 = 0.04, and g2 = 0.1, one gets 

6E 5E' 5E" 

-0.000833975 -0.0008265 -0.0008340 *- 

i.e., an almost perfect agreement, when the In R correction is taken into account. 
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5.3 Double exponential well 



As an example involving smooth potentials, we consider the exponential potential Vi = —gi exp(- 
with gi = I and take ri = 2 for the long-range interaction, and study the changes due to another 
exponential interaction, —A exp(— r/r2) with a much shorter range r2 = 0.02 and a variable 
strength. The results are displayed in Fig. [6} Again, there is a net gain as compared to the ordi- 
nary perturbation theory, and a good agreement with the exact calculation as long as the deviation 
from the unperturbed energy is not too large. 




Xri 



Figure 6: Exponential well supplemented by another exponential of shorter range. Thick 
line: exact, dashed line: SL formula, thin line: perturbation theory. We use here V = 

—gi exp(— r/rl) — A exp(— r/r2) with ri = 2, gi = \ and r2 = 0.02. 



5.4 Harmonic confinement 

The problem is to study how the levels, in particular the ground state, are modified when a 
d = 2 harmonic oscillator is supplemented by a short-range interaction. A recent contribution 
is by Farrell and van Zyl p"4| . They first stressed the property of universality, namely that the 
energy shift does not depend on the details of the short-range potential, but instead is governed 
by the scattering length alone. This is, indeed, a very general property of the exotic atoms, in the 
general sense define in the introduction [j3|. For Vi = r^, the general solution that is regular at 
large distances can be written as 

h{k,r) =exp(-rV2)v^f/(l/2- A;V4,l,r2) , (16) 

in terms of the confluent hypergeometric function U. This expression is simpler, but equivalent to 



the one given in [ 14|. From ([16]), one can calculate explicitly the normalization integral I{E) and 
its derivative. The short-range behavior of h{k, r ) is known and if the ratio of the — to ^/r In r 
coefficients is identified with In a, one recovers the formula given in p4| . Our prescription ( [TT] ) 



corresponds to an approximate, but accurate, solution to the matching equation. For instance, 
using AV2 = —A exp(— r/i?2) as an additional potential, with R2 = 0.02 and A = —80, one 



gets, using the same notation as above and 6E for [ 14|, 



5E 5E 5E' 5E" 

-0.06999 -0.07017 -0.07092 -0.07020 ^ ' 
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Clearly, the main discrepancy comes from reducing this short-range interaction to a zero-range 
ansatz. Once this is accepted, our approximate treatment is nearly exact as compared to the 
precise matching of h(k, r) to the a/t ln(r/a) boundary condition. 



6 Summary 

In this note, we have studied how the energy levels in a wide potential are modified by a short- 
range attraction of increasing strength, focusing on the case of c? = 2 space dimensions, as 
compared to the d = 1 and c/ = 3 situations. 

The energy shifts in a given external potential are well described by the following SL formu- 

r K/a {d = l) 

5E = \ Al/\n{a/R) (d = 2) (18) 
[ Ala {d = 3) 

i.e., a perfect fit is obtained if A^ (and R for = 2) are treated as free parameters. Moreover, 
Aq can be identified with the first non vanishing coefficient of the short-range expansion of the 
radial wave function and is thus proportional to 0(0), the wave function at the origin for the state 
in the external potential alone. The ratio is A^/|0(O)p = 2 (c/ = 1), 2 vr (rf = 2), 47r (rf = 3), 
the unit- sphere area in d dimensions. 

In the d = 2 case which is our main concern, a formula has been derived for R, namely 
InR = —Aq A'{Eo), where A{E) is the coefficient of y/r in the normalized wave function, 
assumed to be real and positive and to match Aq at energy Eq. 

This SL relation becomes more accurate when additional potential V2 becomes more short- 
ranged. In particular, it improves significantly the simple prediction from first order perturbation 
theory in d = 2 and d = 3. 

This study of exotic atoms is intimately linked to the statistical physics of bosons. The 
common tool is the pseudo-potential, which enables one to replace a finite (but short) range 
interaction by a contact interaction. Deriving the pseudo-potential as a function of the scattering 
length for different values of the space dimension d has been extensively discussed. The case of 



c? = 2 is notoriously delicate, see, e.g., p4| - [r7| | for recent contributions. 
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